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ABSTRACT 

We observed 17 optically-selected, radio-quiet high-redshift quasars with the 
Chandra Observatory ACTS, and detected 16 of them. The quasars have red- 
shift between 3.70 and 6.28 and include the highest redshift quasars known. 
When compared to low-redshift quasars observed with ROSAT, these high red- 
shift quasars are significantly more X-ray quiet. We also find that the X-ray 
spectral index of the high redshift objects is flatter than the average at lower 
redshift. These trends confirm the predictions of models where the accretion 
flow is described by a cold, optically-thick accretion disk surrounded by a hot, 
optically thin corona, provided the viscosity parameter a > 0.02. The high red- 
shift quasars have supermassive black holes with masses ~ 10^° Mq, and are 
accreting material at ~0.1 the Eddington limit. We detect 10 X-ray photons 
from the z = 6.28 quasar SDS 1030-1-0524, which may have a Gunn-Peterson 
trough and be near the redshift of reionization of the intergalactic medium. The 
X-ray data place an upper limit on the optical depth of the intergalactic medium 
t{IGM) < 10^, compared to the lower limit from the spectrum of Lya and Ly/3, 
which implies t{IGM) > 20. 

Subject headings: quasars: X-ray - quasars: general - quasars: individual () - 
intergalactic medium - accretion disks 
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1. Introduction 

Rapid growth in the number of quasars discovered with z > 4 has taken place in the last 
few years (Warren et al. 1987; Schneider, Schmidt & Gunn 1989; Kennefick et al. 1995ab; 
Hook, McMahon, Irwin, & Hazard 1996; Storrie-Lombardi, McMahon, Irwin & Hazard 1996; 
Hook & McMahon 1998; Fan et al. 2000, 2001ab; Becker et al. 2001; Schneider et al. 2002). 
We are thus now able to probe the quasar phenomenon in a new regime of high luminosity 
and youth that is important for understanding the structure of quasars and their evolution. 
At the present time, there are about 130 QSOs identified aX z > A (NED 2002), but X-ray 
data for these sources remain sparse, especially for the radio-quiet majority (Bechtold et al. 
1994; Pickering, Impey & Foltz 1994; Henry et al. 1994; Mathur & Elvis 1995; Siebert et 
al. 1996; Fabian et al. 1997; Zickgraf et al. 1997; Siebert & Brinkman 1998; Schneider et 
al. 1998; Wu, Bade & Beckmann 1999; Kaspi et al. 2000; Brandt et al. 2001; Vignah et al. 
2001). The X-ray continuum in active nuclei is thought to arise in a Comptonized wind or 
corona associated with the inner accretion disk around a super-massive black hole (Haardt & 
Maraschi 1993; Czerny & Elvis 1987). Additional components of the X-ray emission are also 
commonly seen in high signal to noise spectra at low redshift, including X-rays associated 
with beamed synchrotron plasma; Compton reflected emission; line emission in Fe Ka; weak 
absorption from warm ionized material; and strong absorption from cold material (see George 
et al. 2000, Reeves & Turner 2000 and references therein). Whether or not there is evolution 
in the structure of quasars with redshift can in principle be learned from X-ray spectroscopy 
of high redshift quasars. 

We report here a survey with the Chandra X-ray Observatory, of optically selected, 
radio-quiet high redshift quasars. This project is part of our long-term program to char- 
acterize the multi-wavelength emission from quasars as a function of redshift (Bechtold et 
al. 1994 ab, Kuhn et al. 2001). The statistical properties of the ensemble population of 
quasars evolve strongly with redshift: the break point, L*, of the luminosity function shifts 
to higher luminosities by a factor of ~50 between z=0 and z=2 (Boyle, Fong, Shanks & 
Peterson 1987). We hope to understand this evolution in terms of a physical model for the 
evolution of individual objects: does this evolution reflect slow changes in a few rare objects, 
or a short-lived phase that most galaxies go through? Are there differences in the spectral 
evolution of radio-loud and radio-quiet quasars? Can these be used to understand the origin 
of the radio emission? The X-ray photons, whose production is closely linked to processes 
in the central engine, should provide direct clues to the physical conditions in the central 
few parsecs, where most of the quasar energy is produced, and will therefore help answer 
these questions. With the recent suggestion that reionization of the intergalactic medium 
happened at 2; ~ 6.3 (Becker et al. 2001), it is possible that 2 ~ 6 is the epoch when the 
first quasars were born. Thus, the quasars targeted in our survey are among the first objects 
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formed in the Universe. 

High redshift quasars are rare, and require pointed observations at known objects for 
study; our expectation is that XMM and Chandra serendipitous surveys will not cover suffi- 
cient area on the sky to find many of them. The increase in sensitivity provided by Chandra 
and XMM however make it possible to increase significantly the number of 2 > 4 quasars 
with X-ray data with only a relatively modest investment of telescope time per object. 

We chose targets from samples of high redshift quasars, based on their optical luminosity. 
The quasars were found in optical surveys, including the APM multicolor survey (Williger et 
al. 1994; Storrie-Lombardi et al. 1994, 1996), the Second Palomar Observatory Sky Survey 
(Smith et al. 1994; Kennefick et al. 1995ab), and the Sloan Digital Sky Survey (Fan et al. 
1999, 2000; Becker et al. 2001). Subsequent observations with FIRST and NVSS confirm 
that are all but one are radio quiet (Stern et al. 2000; Fan et al. 2001b; NED 2002). Since 
these quasars are unusually bright, they have been the subject of other studies, including 
high quahty optical spectrosopy to study the emission lines (Constantin et al. 2001) and 
absorption lines (Storrie-Lombardi et al. 2001; Peroux, Storrie-Lombardi & McMahon 2001). 

We calculated observing times required to detect the quasars in X-rays with 100 photons, 
if the quasar had the average optical to X-ray fiux ratio seen at lower redshift, so that we 
would be assured of a detection (9-10 photons) if the quasar were more X-ray weak. In fact, 
we detected all but one quasar in the sample. Observations of three quasars were taken from 
the public archive, to extend our sample to z = 6.28. 

2. Observations and Analysis 

Seventeen quasars were observed with the Chandra X-ray Observatory and Advanced 
CCD Imaging Spectrometer (ACIS-S, Garmire et al., in preparation; Weisskopf & O'Dell 
1997). All the observations were taken with the quasar on the ACIS-S3 CCD, and reduced 
with the standard pipe-line reduction software, and CIAO (2002; version 2.2). Initially, we 
pointed off-axis to mitigate pile-up in anticipation of high count rates, but later moved the 
quasar position on axis. Pile-up is negligible for all observations in our survey. 

The observations are summarized in Table 1. We list the exposure times, observation 
times and net counts detected. For every object except PSS 1435+3057, we detected a 
significant X-ray source within 1 arcsec of the optical position, so there is no doubt to the 
identification. 

For PSS 1435-1-3057, we give a 3a upper limit to the X-ray photons at the optical 
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position. There is a weak source (11 photons) located 3.8 arcsec from the optical position. 
This source is beyond the range of aspect errors (Aldcroft 2002), but may be associated with 
an extended structure from the quasar (c.f. Schwartz 2002). 

In Table 2, we list X-ray flux rates and other X-ray parameters. Although the number 
of photons detected is small, we were able to derive meaningful spectral fits. We used Sherpa 
(CIAO 2002, version 2.2) to fit power law parameters to the counts, including absorption by 
the Milky Way column density of neutral hydrogen fixed at the value inferred from 21-cm 
maps (GOLDEN 2002, which is based on Dickey & Lockman 1990). We fit a function of the 
standard form 

where is the energy index and fo is the normalization at 1 keV, with units of photons keV~^ 
cm~^ s~^. We restricted our fits to the range 0.3-6.5 keV, where the calibration is reliable at 
the present time, and the Chandra background is the lowest. The source extraction region 
for each quasar was a circle with 10 pixels (4.92 arcsec) radius. There were no background 
flares during the observations. 

The choice of cosmology is important when converting fluxes to luminosities. The 
absolute B magnitude, Mb, is related to apparent magnitude, niB, by 

MB = mB + 5-5 log D (2) 

where 

D = ^A (3) 

where c is the speed of light and Hq is the Hubble constant. 

In general (Peacock 1999, his equation 3.39), for a flat Universe, 

A = {l + z) [ {nA + nM{l + zf)-^^^ dz (4) 
Jo 

For go > and fiA=0, we have 

A = z{l + '-^^^ r) (5) 

^ (l + 2go^)0-5 + (l + go^)^ 
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which simphfies to 

A = z{l+'-) (6) 
for go = (see also Carroll, Press & Turner 1992). 

The quasar luminosity function is usually reported assuming go = V^? Ho= 50 km 
Mpc~^, A=0, while observers often adopt go = and Hq = 70 km Mpc~^ and A=0. 
The "best" fit to cosmic microwave background data, Type la SNe light curves, and large 
scale structure models suggest that Qa = 0.7, Qm=0.3, and Ho= 70 km s~^ Mpc~^. We 
refer to this choice of cosmological parameters as "ACDM" in calculations below. Note that 
the integral of equation (4) must be evaluated numerically. Below we report luminosities for 
each of these three cosmologies. 

3. The Evolution of ctox for Radio-Quiet Quasars 

Following Zamorani et al. (1981) we compute a ratio of X-ray to optical flux, aox, where 

ttox = -logifx/ fopt) /logiu^/uopt) (7) 
and log Ux = 17.6845 for rest-frame 2 keV and log Uopt = 15.0791 for rest frame 2500 A. 

In Table 2 we give the results. For all objects, the continuum flux at 1450 A is available 
from the literature (see Table 1), usually in the form of apparent AB magnitude at 1450 
A in the rest frame of the quasar, m^^, where 

mAB = -2.5 logU - 48.57 (8) 

so that fi, has units of ergs cm~'^ Hz~^. We compute fi, at 2500 A, which has units 
of ergs s~^ cm'"^ Hz~^, assuming the optical continuum has a power law spectral energy 
distribution with frequency spectral index a = -0.3. We compute fx at 2 keV in the rest 
frame of the quasar from the measured flux at 1 keV in our observed frame, assuming a 
power law with energy index rx=2.2. 

Figure 1 shows aox versus redshift, z, and Figure 2 shows aox versus absolute B mag- 
nitude, Mb- For clarity, we plot the errors on aox only for the quasars with z > 5. For the 
rest, although the number of X-ray photons detected is small and the uncertainties on fx 
correspondingly large, the division by log{vx/ ^opt) = 2.60 makes the errors only ~ 0.2. 
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For a low redshift comparison sample, we use the compilation of ROSAT all-sky survey 
and pointed observations for radio-quiet quasars (Yuan et al. 1998) and radio-loud quasars 
and blazars (Brinkman et al. 1997), supplemented with observations of other high redshift 
quasars from the literature (Schneider et al. 1998; Kaspi, Brandt & Schneider 2000; Brandt 
et al. 2001; Vignah et al. 2001). Since Yuan et al. (1998) and Brinkman et al. (1997) 
did not list Oox, we computed it from the listed unbsorbed X-ray flux density (0.1-2.4keV), 
assuming that the spectrum is a power law, with energy index with the best fit listed 
by Yuan et al. (1998). For objects with no F measured, we assume F^. = 2.58 for z < 0.5, 
F^ = 2.46 for 0.5 < z < 1.0, F^ = 2.35 for 1.0 < ^ < 2.0 and F^ = 2.2 for z > 2.0 (Yuan et 
al. 1998). For the optical flux, we used the absolute B-magnitude for each object listed in 
Veron-Cetty & Veron (2000) which includes a K-correction. We extrapolated to 2500 A in 
the rest frame of the quasar assuming a power law with a = -0.3. The data points plotted 
in Figure 1 and Figure 2 are available in the electronic version of this paper, and on our 
website. 

The z > 4 radio-quiet quasars are clearly more X-ray quiet than their low-redshift 
counterparts, even when their extreme luminosity is taken into account. Previous studies 
(e.g. Avni, Worrall & Morgan 1995; Brandt et al. 2001; and references therein) found that 

depends mostly on optical luminosity, although the correlation of luminosity and redshift 
in the observed samples made it difficult to sort out whether Oox depended mostly on optical 
luminosity or redshift (e.g. Bechtold et al. 1994; Pickering, Impey & Foltz 1994). The 
Chandra sample also suffers from a strong correlation between redshift and luminosity. In 
Figure 1 the three Chandra quasars at z ~ 3 are strikingly offset from the ROSAT quasars to 
larger aox, while in figure 2 they show no offset. This is likely due to the newer surveys from 
which they are taken which cover larger solid angle than those used to select the ROSAT 
high z quasars. As a result they find systematically more luminous quasars. 

To quantify the result, we computed the generalized Kendall's tau, including lower 
limits on aox for quasars not detected at high redshift (Avni 1976; Avni, Soltan, Tananbaum 
& Zamorani 1980; Feigelson & Nelson 1985; Avni & Tananbaum 1986; Isobe, Feigelson & 
Nelson 1986; Akritas & Siebert 1996). We used the IRAF program bhkmethod to search 
for correlations between (1) redshift and aox and (2) absolute B magnitude and aox- 

We find that aox is anticorrelated with redshift, with Kendall's tau = -0.38, and Z-value 
=8.4, so that a correlation is present at 8Aa significance. On the other hand, the probability 
is 0.05 that aox is not correlated with absolute B-magnitude (Z-value=1.9). Thus, we find 
that aox depends primarily on redshift, and weakly on luminosity. 

Does this conclusion depend upon our assumptions for the X-ray spectral index, F^^? 
In figure 4, we plot F^; versus redshift and versus luminosity, for the Yuan et al. (1998) 
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ROSAT sample and our high redshift sample. The high redshift quasars if anything are 
flatter (smaller F^.) than their low redshift counterparts. If we had adopted a mean of 
1.5 instead of 2.2, we would have increased by ~0.3, making the difference between low- 
and high- redshift greater. Thus the assumed spectral indices cannot account for the result. 

Figures 1 and 2 confirm the result seen previously that radio-loud quasars are more 
X-ray loud than radio quiet quasars, at all redshifts and luminosities. 

Although we have emphasized how X-ray quiet the high redshift quasars are given their 
optical luminosity, they are still prodigious emitters of X-ray photons. In Table 3, we list 
the X-ray luminosity (erg s~^), in the 2-10 keV band, computed from the fits in Table 2. 
The quasars have X-ray luminosities of 10^^ — 10^^ erg s~^. 

4. Comparison with Other Studies 

After submission of this paper, three papers were submitted to the Astrophysical Jour- 
nal Letters about the Chandra observations of SDS 0836-1-0054, SDS 1306-1-0356 and SDS 
1030-1-0524. These quasars were observed with director's discretionary time and were placed 
in the public archive immediately. Brandt et al. (2002) and Mathur, Wilkes & Ghosh (2002) 
report fluxes for the three quasars and conclude, contrary to the results reported here, that 
the optical-to-X-ray flux ratios of these three highest redshift quasars are not significantly 
different from those of low redshift quasars. The different conclusion is a result of the dif- 
ferent comparison samples used, which in both papers were smaller than the one presented 
here. Schwartz (2002) reported not only the core fluxes for the three quasars, but also argued 
that an X-ray source 23" from SDS 1306-1-0356 with no optical counterpart on the Palomar 
Sky Survey is in fact associated with the quasar. He postulates that the X-rays are Compton 
scattered cosmic microwave background photons from a jet structure, similar to those seen at 
lower redshift (Tavecchio et al. 2000; Celotti, Ghisellini & Chiaberge 2001; Siemiginowska et 
al. 2002). Future observations can confirm this intriguing result. In our sample, as discussed 
above, we see a similar source 3.8" from the optical position of PSS 1435-1-3057. 

For the core X-ray fluxes, we measure slightly different values than those reported in the 
three other papers. We measure 24, 10 and 19 photons for SDS 0836-1-0054, SDS 1306-1-0356 
and SDS 1030-1-0524, respectively. Schwartz (2002) measures 21, 6, 16 photons; Mathur, 
Wilkes & Ghosh (2002) measure 21, 6, 18 photons; and Brandt et al. (2002) measure 21, 6 
and 17 photons. The difference lies in two factors: the energy range used (we use 0.3 - 6.5 
keV, whereas the others use 0.5 - 7 keV or 0.5 - 8 keV) and the difference in circle radius 
used for source extraction (we use 4.9 ", whereas they use 1.2 - 2.9 "). We chose the 0.3 - 6.5 
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keV range because it has the lowest background, and the larger circle extraction to include 
the point spread function at all energies. The result is a greater number of net photons for 
all sources, and a more reliable measure of the source properties. 

5. Mean X-ray Spectrum 

Figure 3 suggests that the X-ray spectral index flattens at high redshifts. This could 
be caused by two effects. First, as the redshift increases, the observed X-ray band samples 
greater and greater energies. Thus, the flattening may be caused by intrinsic flattening of 
the power law at high energies. Second, the quasars in our survey may have flatter intrinsic 
spectra at all energies, due to their high luminosities, or high redshift. 

We applied Kendall's generalized tau to investigate the dependance of F^; on redshift 
and optical luminosity. We did not include objects for which a power law index had not 
been measured. We found that F^, is anti-correlated with redshift (Kendall's tau=-0.4319, 
Z-value = 8.46) and positively correlated with luminosity (Kendall's tau=0.2330, Z-value = 
4.57). The data suggests that F^. depends strongly on both redshift and luminosity. 

We note that at the redshifts of the quasars in our sample, the usable Chandra energy 
range, 0.3-6.5 keV, corresponds to relatively hard X-rays, 1.5 keV to 32 keV. Therefore, we 
expect absorption by intervening damped Lya systems and other intervening absorbers to 
be neglible, unless the absorbers are at low redshift. Since we chose the sample objects to 
avoid known broad-absorption line quasars (see Green et al. 2001), intrinsic absorption is 
also probably negligible. We discuss the special case of SDS 1030+0524 at z = 6.28 below. 

6. Implications for Mbh and rii 

A successful model for explaining the optical/ultraviolet continuum of quasars (the Big 
Blue Bump) involves accretion onto a supermassive black hole through an optically thick, 
physically thin disk (Shields 1978; Malkan 1983; Bechtold et al. 1987; Czerny & Elvis 1987; 
Sun & Malkan 1989). The emitted spectrum is then the integral of Planck spectra of different 
temperatures, resulting in a fiat continuum, oc through the optical-UV. There is a near 
exponential falloff blueward of a cutoff energy, Eco, since there is a maximum temperature for 
the accreting material, nearest the black hole. For the supermassive black holes in quasars 
accreting near the Eddington limit, the cutoff is Eco ~ 10-100 eV. 

The spectral shape and luminosity of the accretion disk can be predicted and depends 
on five parameters: the black hole mass (Mbh)', the mass accretion rate m; the total angular 
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momentum; the viscosity in the disk, a; and the inchnation at which we observe the disk. 
If we make the usual assumptions that the black hole is maximally spinning (Kerr black 
hole; see Elvis, Risalti & Zamorani 2002), that we are observing the disk face-on, and that 
the viscosity a=0.1, then the spectral energy distribution and luminosity depends on two 
parameters, Mbh and m. 

Malkan (1990) showed that for a Kerr black hole, an analytical expression for the emitted 
spectrum is 

A{—Y/'exp{-—) (9) 

^co ^co 

where A and are the normalization and cutoff frequency. The cut-off energy can then be 
written (Wandel 2000) as 

hu,o = {6eV)m^/^M-^^^ = {20eV) lH^M'^^^ (10) 

where Mg is the black hole mass, in units of 10^ solar masses. 

Observationally, is it also possible to relate Mbh to the width of the broad emission 
lines, assuming the BLR gas is in Keplerian motion, and calibrating the size (and hence 
ionization parameter) of the gas by reverberation mapping of low-redshift AGN (Peterson & 
Wandel 1999; Wandel, Peterson & Malkan 1999). 

Thus, with two observations - the width of H/? and the continuum luminosity - we 
can solve for the two free parameters of the model, Mbh and rh (Wandel & Petrosian 1988; 
Wandel & Boiler 1998; Wandel 1999). In Figures 4 and 5 we plot the result. For comparison, 
we also plot the results for narrow line Seyfert I's (Crenshaw 1986; Stirpe 1990; Boiler, 
Brandt & Fink 1996; Brandt, Mathur & Elvis 1997), the PG quasars (Boroson & Green 
1992; Miller et al. 1992), and the LBQS quasars (Forster et al. 2001). For high redshift 
quasars, B.j3 is in the near-IR, and only the most luminous quasars have measurements. We 
use H/3 whenever possible (Hill, Thompson & Elston 1993; Rokaki, Boisson, & Collin- Souffr in 
1992; Nishihara et al. 1997; Mcintosh et al. 1999). For the rest, we use the FWHM of the 
CIV emission line, and assume the relation given by Corbin (1991), 

FWHM{Hf3) = 1.35 FWHM{CIV) - 1391 km s'^ (11) 

We use C IV for intermediate redshift quasars from the LBQS (Forster et al. 2001), and 
observations of 2 ~ 4 quasars from Constantin et al. (2002). To estimate the luminosity, we 
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used the optical luminosities computed from the magnitudes listed in Table 1, or from the 
Veron-Cetty & Veron (2000) catalog. For all quasars, we apply a bolometric correction of 
10. 

For the ACDM cosmology (Figure 5) we see that the NLSyl's have relatively low mass 
black-holes, PG quasars are more massive, and the high redshift quasars have very massive 
black holes. The typical black hole mass for the z = A — 6 quasars is ~ 10^° Mq, with most 
objects accreting between 10 ~^ and 0.8 of the Eddington mass accretion rate. 

7. Accretion disk and X-ray emitting corona 

We now use the X-ray measurements to further investigate the nature of the black hole 
and accretion in quasars. Janiuk & Czerny (2000) have presented calculations of the X-ray 
spectrum expected from a hot corona associated with an accretion disk. The model predicts 
the fraction of energy dissipated in the corona as a function of the disk radius. The coronal 
dissipation is assumed to be proportional to the gas pressure, and the pressure at the base of 
the corona is determined by the condition for the disk/corona transition. This requirement 
is consistent with the evaporation/condensation equilibrium condition (Rozanska & Czerny 
2000) and effectively means that bremsstrahlung and Compton cooling at the base of the 
corona are comparable. The spectrum from the disk and the corona are computed locally, 
taking into account the Comptonization in the hot coronal plasma. 

The observed model spectrum is integrated over the disk surface. The spectrum is 
parameterized by the black hole mass, accretion rate and coronal viscosity, a. Here we used a 
code which includes tabulated amplification factors for Comptonization in the corona derived 
from Monte Carlo simulations (Janiuk, Czerny & Zycki 2000), instead of the analytical 
approach used by Janiuk & Czerny (2000). In contrast to the calculations described in 
the previous section, we assume a Schwarzschild, not Kerr black hole. We calculate model 
spectra assuming the disk is face-on. Based on low luminosity accretion models (ADAFs, 
Kurpiewski & Jaroszynski 2000, 1999) we expect that the X-ray spectra will be harder for 
the maximally rotating black hole. Kerr geometry should be considered in the future models 
to quatify the effects. 

We defer a detailed comparison of the new data with the models to a future paper, 
and here present only the general trends. In Figure 6 we plot representative spectral energy 
distributions in order to illustrate the dependence on parameters. We plot models with the 
black hole mass fixed at 10^*^ M© and different viscosities and accretion rates. We see that 
the X-ray spectral index is only a weak function of the parameters, but that the optical/UV 
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luminosity and aox are strong functions of the viscosity, a, and accretion rate. We use the 
standard a viscosity prescription (Shakura & Sunyaev 1973) which describes the efficiency 
of angular momentum transfer in the disk. 

Figure 7 shows how aox and depend on parameters. We plot the results for models 
with black hole mass 10'' and lO^'' and accretion rates of 0.01 and 0.8 Eddington, which 
bracket the values derived for the luminous quasars in Figures 4 and 5. 

In Figure 7 we see the following trends: 

1. The large observed scatter in aox is predicted naturally by the models, for a plausible 
range of accretion rate. 

2. If the luminous high redshift quasars have on average more massive black holes than 
lower redshift quasars, then they are predicted to have larger aox- That is, the quasars in 
our sample are predicted to be more X-ray quiet than the PG or LBQS quasars, in agreement 
with the observations. 

3. For the massive black holes in the quasars in our sample, aox ~ 2 implies that the viscosity 
parameter a > 0.02. 

4. The predicted X-ray spectral index, T^, is between 1.2 and 1.6, similar to what we measure 
for the high redshift sample. The low redshift quasars with steeper F^, must have an extra 
soft component, or more complex spectra than the simple power law fits. 

The models predict the spectrum throughout the ionizing ultraviolet and soft X-rays, 
and therefore can be used to predict emission line properties, particularly the fiux and 
equivalent width of C IV. The implications for the broad line region will be discussed further 
in a future paper. 

8. Detection of SDS 1030+0524 at z = 6.28 

SDS 1030+0524 is the highest redshift quasar discovered to date. Becker et al. (2001) 
present the ultraviolet spectrum which shows very strong absorption just blueward of the 
Lya emission line. They convincingly argue that the absorption is far stronger than what 
is predicted from a simple extrapolation of the Lya forest from lower redshift, so that they 
have detected a Gunn-Peterson trough. They conclude that the intergalactic medium (IGM) 
was in the process of reionization at 2; ~ 6. They derive a lower limit on the optical 
depth of neutral hydrogen in the intergalactic medium, from analysis of Lya and Ly/5, of 
r{IGM) > 20. 
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With Chandra, we detect 10 photons for SDS 1030+0524 in the observed 0.3-6.5 keV 
band, which is 2.2-47.0 keV in the quasar's rest frame. In fact all 10 photons have E < 2.5keV 
observed, ot E < 18 keV. The limits on IGM absorption are weak. We used Sherpa to fit 
the 10 photons with a model which had the Galactic absorption fixed, a power law with 
fixed r^,. = 2.2, and absorption at z = 6.28. The column of the redshifted absorption and 
normalization of the flux were free parameters. Significant absorption was not detected, and 
we can place a 3a upper limit to the absorbing column of 7.6x10^^ atoms cm~^ if the IGM 
has solar abundance, and 5.3x10^^ atoms cm~^ if the IGM has hydrogen and helium only, 
at the solar ratio. In either case, the the Chandra data imply that at the Lyman limit, 
t{IGM) < 10^ at redshift ~ 6. 



9. Discussion 

We have studied the evolution of quasars from 2; = 6 to the present day, putting together 
what is known about the optical continuum luminosities, broad emission line widths, and 
X-ray flux and spectrum. The data are consistent with a model where the optical-ultraviolet 
continuum arises in a cold, thin, optically thick accretion disk, and the X-rays are produced 
at high redshift by a Comptonized corona. At low redshift, extra components are likely 
making the X-ray spectrum more complex than the simple power law fits available so far for 
most objects. 

The high redshift quasars have systematically more massive central black holes than 
their low-redshift counterparts, and are accreting at high rates, several tenths of the Edding- 
ton limit. Their observed relatively weak X-ray fluxes are a natural outcome of the accretion 
disk-coronal models, with no cold absorption necessary to suppress the X-ray emission (c.f. 
Brandt et al. 2001; Mathur 2001). 

The Chandra observations require large values of the viscosity parameter, a > 0.02 for 
the high redshift quasars. Numerical simulations show that the turbulent part of the a pa- 
rameter is negligible in comparison to the magnetic term arising from magnetorotational in- 
stability and a should be within 0.001-0.1 (Balbus & Hawley 1998, Armitage 1998, Armitage 
et al 2001). Observational constraints based on AGN variability are limited (Siemiginowska 
& Czerny 1989), but in general they are in agreement with the theoretical predictions, al- 
though for high disk luminosities a exceeds 0.1 for the PC sample of quasars (Starling et al 
2002). 

We conclude that even at the highest redshifts probed to date, quasars were producing 
prodigious ultraviolet and soft-X-ray photons, which no doubt had interesting effects on the 
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intergalactic medium and formation of galaxies at the earliest times. 
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Table 1. Summary of Observations 





Quasar 




AB" 


ObsID 


Exp. 


RA,Dec 


^^(Gal) 


Net 


Refn. 










Date 


(sec) 


J2000 


Counts'^ 




PSS 


0059+0003 


4.178 


19.45 


2179 
18Sep2001 


2682 


00:59:22.80 
+00:03:01.0 


3.20 


12 


1,2 


BRI 


0103+0032 


4.437 


18.84 


2180 
18Sep2001 


3709 


01:06:19.20 
+00:48:22.0 


3.10 


26 


1,2,3 


SDS 


0150+0041 


3.67 


18.35 


2181 
3lAug2001 


3238 


01:50:48.80 
+00:41:26.0 


2.79 


7 


8 


BRI 


0241-0146 


4.053 


18.45 


875 
llMar2000 


7365 


02:44:01.90 
-01:34:03.0 


3.73 


17 


4 


PSS 


0248+1802 


4.43 


18.24 


876 
27Decl999 


1731 


02:48:54.30 
+18:02:50.0 


9.72 


19 


1,2 


BRI 


0401-1711 


4.236 


18.84 


2182 
3Aug2001 


3841 


04:03:56.60 
-17:03:24.0 


2.34 


15 


4 


SDS 


0836+0054"' 


5.82 


18.8 


3359 
29Jan2002 


5687 


08:36:43.85 
+00:54:53.3 


4.15 


24 


9 


SDS 


1030+0524 


6.28 


19.7 


3357 
29Jan2002 


7955 


10:30:27.10 
+05:24:55.1 


3.09 


10 


9 


BRI 


1033-0327 


4.509 


18.84 


877 
26Jan2000 


3447 


10:36:23.70 
-03:43:20.0 


4.85 


16 


5,6 


PSS 


1057+4555 


4.10 


17.53 


878 
14Jun2000 


2808 


10:57:56.40 
+45:55:52.0 


1.17 


34 


2 


SDS 


1204-0021 


5.10 


19.05 


2183 
2Dec2000 


1570 


12:04:41.70 
-00:21:49.0 


2.14 


26 


7 


SDS 


1306+0356 


5.99 


19.6 


3358 
29Jan2002 


8156 


13:06:08.26 
+03:56:26.3 


2.08 


19 


9 


PSS 


1317+3531 


4.36 


19.55 


879 
14Jun2000 


2788 


13:17:43.20 
+35:31:31.0 


0.99 


9 


2 


PSS 


1435+3057 


4.35 


19.12 


880 
21May2000 


2811 


14:35:23.50 
+30:57:23.0 


1.21 


< 9 


1,2 


PSS 


1443+2724 


4.42 


19.23 


881 
12Jun2000 


2170 


14:43:31.20 
+30:57:23.0 


2.33 


10 


1,2 
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Table 1 — Continued 



Quasar 






ObsID 


Exp. 


RA,Dec 


^^(Gal) 


Net 


Refn. 








Date 


(sec) 


J2000 


Counts^ 




SDS 1621-0042 


3.70 


17.41 


2184 


1570 


16:21:16.90 


7.29 


26 


7 








5Sep2001 




-00:42:51.1 








BRI 2212-1626 


3.99 


18.65 


2185 


3222 


22:15:27.20 


2.64 


15 


4 








16Dec2001 




-16:11:33.0 









'^AB magnitude at 1450A in quasar restframe, AB = -2.5 log fopt - 48.57, where [fopt] = 
erg cm~^ Hz~^. 

^Galactic H I column, in units of 10^° atoms cm~^, from GOLDEN. 
'^Net number of photons in observed energy range 0.3 - 6.5 keV 
<^SDS 0836+0054 is radio-loud - see Fan et al. (2001b). 

References. — (1) Kennefick, J. D. et al. (1995); (2)Kennefick, J. D. et al. (1995); 
(3)Smith,J. D., et al. (1994); (4)Storrie-Lombardi,L. J., et al. (1996); (5)Storrie-Lombardi,L. 
J., et al. (1994); (6)Williger, G. M., et al. (1994); (7)Fan, X., et al. (2000); (8)Fan, X., et 
al. (1999); (9) Becker et al. (2001). 
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Table 2. X-ray Results 



QSO 






pa 




J O 


log f2keV 


log /2500 




a'' 

"ox 


log Fl 


rbbUUoy+UUUo 


A 1 TO 


\j. 


_i_n 

.75±0 


A 1 

.41 


O O O _l_ 1 O [X 

3.38±1.25 


-31 


150 


-27 


137 


1, 


t-/|r|+0.276 
•^^U_o.236 


1 1 T 

-13.17 


"DT3Tm no 1 nnoo 

BKi(Jl(Jo+(J(Joz 


A AO'T 

4.4d7 


1, 


TO _i_n 


.30 


T OO A O 


-30 


803 


-26 


893 


1, 


t-rii +0.245 
•^Ui_o.221 


10 CO 


oDbUioU+UU41 


O.O { 


0, 


oo 1 n 

.oo±U 


.4 i 


1.1 ±U.44 


-31, 


.668 


-26 


697 


1, 


Qi 1 +0.259 
•y-"- -'--0.219 


1 0^ 


l3KiUz41-U14b 


4.U0O 


i, 


oo 1 n 

.oo±U 


.36 


z.oo±U.oy 


-31, 


.287 


-26 


737 


i. 


74Q+0.259 
• '^^-0.219 


1 TO 


T~)CCno/io 1 1 ono 

raoU248+l8U/ 


A AO 

4.4o 


1, 


c\o _i_n 

.93±0 


.37 


1 A TO _l_0 yl O 

14.78±3.48 


-30, 


.488 


-26 


653 


1, 


/I7C-+0.252 
'-'-0.226 


1 00 

-13.38 


"DTDTri/im 1 '71 1 


4.23d 


1, 


OO _i_n 

.28±0 


.37 


o ^n_i_i nn 

3.60±1.00 


-31, 


.121 


-26, 


,893 


1, 


(-r,f-+0.255 
•"^^-0.226 


10 /I 

-13.54 


bUD(Jooo+UU54 


5.82 


1, 


/I /I 1 r\ 

.44±0 


on 

.30 


/I on_i_n no 

4.30±0.93 


-30, 


.906 


-26, 


,877 


1, 


r/(Q+0.290 

•^^^-0.260 


10 C T 

-13.57 


SDS1030+0524 


6.28 


1, 


.36±0 


.46 


1.21±0.58 


-31, 


.422 


-27, 


,237 


1, 


•609+°io'5 

— U.oUO 


-13.07 


BR11033-0327 


4.509 


2, 


.39±0 


.42 


5.01±1.29 


-30, 


.951 


-26. 


,964 


1, 


t-f.Q+0.260 
•^""-"-0.232 


-14.14 


PSS1057+4555 


4.10 


1, 


.80±0 


.26 


11.47±1.97 


-30, 


.631 


-26. 


,369 


1, 


(-00+0.228 
•"'-'°-0.207 


-13.39 


SDS1204-0021 


5.10 


1, 


.50±0 


.38 


2.30±0.61 


-31, 


.236 


-26. 


,977 


1, 


(-07+0.280 
•"'J ' -0.250 


-13.89 


SDS1306-f0356 


5.99 


1, 


.32±0 


.33 


2.15±0.52 


-31, 


.194 


-27. 


,197 


1, 


cr 07+0.300 
•^'J ' -0.269 


-13.79 


PSS1317+3531 


4.36 


2, 


.36±0 


.57 


2.83±1.01 


-31, 


.213 


-27. 


,177 


1, 


crcr9+0.279 
•"-'"-'^-0.240 


-14.36 


PSS1435+3057 


4.35 








<3.0^' 


<-31.189 


-27, 


,005 




>1.60 




PSS1443+2724 


4.42 


1, 


.50±0 


.47 


4.47±1.46 


-31, 


.009 


-27, 


,049 


1, 


[-99+0.273 
•'-'^'^-0.238 


-13.60 


SDS1621-0042 


3.70 


1, 


.39±0 


.30 


18.34±4.02 


-30, 


.469 


-26. 


,321 


1, 


p-qc:+0.222 
•'-'^'-'-0.200 


-12.91 


BRI2212-1626 


3.99 


1, 


.27±0, 


.37 


4.32±1.21 


-30, 


.724 


-26, 


,817 


1, 


^04+0.247 
•"<J^-0.218 


-13.45 



= power law energy index, A{E) = fo[E/lkeV] Absorption fixed at Galactic NH 

^fo = normalization of power law, in units of 10"^ photons keV~^ crn"^ at 1 keV. 
Note that f^ (ergs cm'^ sec"^ }iz-^)=6.63xl0~^'^ foE{keV). 

^hkev = extrapolated flux at E=2 keV in the quasar rest frame, erg cm~^ s^^ Hz~^, 
assuming r2;=2.2. 

'^/25oo = extrapolated flux at 2500A, calculated from measured flux at 1450A, assuming 
power law with a=-0.3; erg cm~^ s~^ Hz~^. 

'^O-ox = ratio of X-ray to optical flux, see equation (7) in text. Errors include uncertainties 
from tabulated uncertainties in fo and P^,. only. 

^Integrated observed flux within 2-10 keV, ergs cm~^ s~^. 
^Three sigma upper limit, assuming P^; = 2.2. 
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Table 3. X-ray Luminosities 



Quasar 




LogLx 


1 r\n T 

LOgLx 


lr\n T °' 






go — U 


„ 1 /9C 


ivC- U IVl 


PSS 0059+0003 


4.178 


46.39 


45.70 


46.06 


BRI 0103+0032 


4.437 


46.07 


45.35 


45.71 


SDS 0150+0041 


3.67 


46.11 


45.48 


45.83 


BRI 0241-0146 


4.053 


45.74 


45.06 


45.42 


PSS 0248+1802 


4.43 


46.27 


45.54 


45.91 


BRI 0401-1711 


4.236 


46.04 


45.34 


45.70 


SDS 0836+0054 


5.82 


46.48 


45.62 


46.00 


SDS 1030+0524 


6.28 


47.10 


46.19 


46.58 


BRI 1033-0327 


4.509 


45.53 


44.80 


45.17 


PSS 1057+4555 


4.10 


46.14 


45.46 


45.82 


SDS 1204-0021 


5.10 


45.97 


45.17 


45.54 


SDS 1306+0356 


5.99 


46.30 


45.43 


45.81 


PSS 1317+3531 


4.36 


45.26 


44.55 


44.91 


PSS 1443+2724 


4.42 


46.04 


45.32 


45.69 


SDS 1621-0042 


3.70 


46.48 


45.84 


46.20 


BRI 2212-1626 


3.99 


46.04 


45.37 


45.73 



^X-ray luminosity, in observed 2-10 keV band, in ergs 



''qo=0, Ho=70 km s"^ Mpc-\ A = 0. 
'=qo=l/2, Ho=70 km s'^ Mpc-\ A = 
^^A = 0.7, VLm = 0.3, Ho=70 km s~^ Mpc" 
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Fig. 1. — X-ray to optical flux ratio versus redshift, z, for radio-quiet (top) and radio- loud 
(bottom) quasars. X's are quasars detected by ROSAT, from Yuan et al. (1998) for radio 
quiet objects and Brinkman et al. (1997) for radio loud quasars; Open circles are high 
redshift quasars observed with Chandra (see text). 
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Fig. 2. — X-ray to optical flux ratio versus absolute B-magnitude, M^, computed assuming 
ACDM cosmology. X's are radio quiet quasars from Yuan et al. (1998) and radio-loud 
quasars from Brinkman et al. (1997). Open circles are high redshift quasars observed with 
Chandra (see text). 
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6 



Radio Quiet Quasars 




Fig. 3. — X-ray spectral index, F^, versus versus redshift, z. X's are radio quiet quasars 
from Yuan et al. (1998); Open circles are high redshift quasars (see text). For quasars with 
too few counts detected to derive F^, we plot F2-=0. 
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n I \ I I I \ I I I \ I I I ^ I I I ^ I I r 

Black Hole Parameters 




-20 -22 -24 -26 -28 -30 -32 

Mg (qo=l/2, Ho=50 km s-i Mpc-i) 



Fig. 4. — Black hole parameters. Plotted is the UV cut-off energy, Eco in eV, versus absolute 
B magnitude, M^, for go = 0.5, A=0, and Hq = 50 km s~^ Mpc~^. Small open circles are 
Narrow Line Syl's, triangles are PG quasars, open squares are LBQS quasars, and large 
open circles are high redshift quasars. Dotted lines show locus of constant black hole mass, 
for Mbh = 10^ 10^ 10^ 10^ 10^°, and 10^^ M©. Long dashed lines are locus of constant 
mass accretion rate, in units of LO, 0.1, 0.01 and 0.001 the Eddington rate. 
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n 1 ^ 1 1 1 ^ 1 1 1 ^ 1 1 1 ^ 1 1 1 ^ 1 1 r 

Black Hole Parameters 




-20 -22 -24 -26 -28 -30 -32 

Mg (ACDM) 

Fig. 5. — Same as figure 4, except Mb was computed assuming ACDM cosmology. 
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\ \ 
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Fig. 6. — Representative spectral energy distributions for two-phase accretion disk models of 
Janiuk & Czerny (2000). We assume a 10^'' Mq black hole, and viscosity parameter a=0.01 
and 0.5, for mass accretion rate = 0.8 and 0.01 times the Eddington limit. 
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Fig. 7. — Predicted dependence of aox (top panel) and X-ray power law spectral index 
(bottom panel) on viscosity, a, for two-phase accretion disk models of Janiuk & Czerny 
(2000). Dotted lines are for a supermassive black hole with 10^" Mq, dashed lines have 10^ 
Mq black hole. Points are plotted for accretion rates of 0.01, 0.1, 0.3, 0.5, 0.8 times the 
Eddington limit. 



